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ABSTRACT 



7 We present an analysis of the gamma-ray measurements by the Large Area 
Telescope (LAT) onboard the Fermi Gamma-ray Space Telescope in the region of 
the supernova remnant (SNR) Cygnus Loop (G74. 0—8.5). We detect significant 
gamma-ray emission associated with the SNR in the energy band 0.2-100 GeV. 
The gamma-ray spectrum shows a break in the range 2-3 GeV. The gamma-ray 
luminosity is ~ 1 x lO^^erg s~^ between 1-100 GeV, much lower than those of 
other GeV-emitting SNRs. The morphology is best represented by a ring shape, 
with inner/outer radii 0°.7 ± 0°.l and 1°.6 ± 0°.l. Given the association among 
X-ray rims. Ha filaments and gamma-ray emission, we argue that gamma rays 
originate in interactions between particles accelerated in the SNR and interstel- 
lar gas or radiation fields adjacent to the shock regions. The decay of neutral 
pions produced in nucleon-nucleon interactions between accelerated hadrons and 
interstellar gas provides a reasonable explanation for the gamma-ray spectrum. 

8 Subject headings: cosmic rays — acceleration of particles — ISM: individual objects 

9 (the Cygnus Loop) — ISM: supernova remnants — gamma rays: ISM 
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Introduction 



Diffusive acceleration by supernova shock waves can accelerate particles to very high 



12 energies (e.g. 



Blandford fc Eichler 



19871 ). Gamma-ray observations are a useful probe of 
these mechanisms complementary to other wavebands. So far, observations by the Large 
Area Telescope (LAT) on board the Fermi Gamma-ray Space Telescope have demonstrated 



that bright gamma-ray sources coincident with midd 



16 interacting with dense molecular clouds ( lAbdo et al. 



e-aged supern ova remnants (SNRs) 



2009 



2010a 



exhibit steep 



17 gamma-ray spectra above a few GeV. A possible conventional explanation for these spectral 



IB properties is that the energy distribution of cosmic 



19 their diffusive transport (e.g 



2010 



Ohira et al. 



Aharonian fc Atoyan 



rays 



1996 



CR) is greatly influenced by 



Gabici et al. 



2009 



Torres et al. 



201ll ). On the other hand, these features can also be explained by 



reacceleration of pre-existing cosmic rays at a cloud shock and subsequent adiabatic 
compression where strong ion-neutral collisions accompanying Alfven wave evanescen ce 



23 lead to a steepening of the spectrum of accelerated particles (jUchiyama et al. 



24 Furth ermore, using three-dimensional magnetohydrodynamic simulations. 



2010h 



Inoue et al 



25 ( I2OIOI ) show that the interaction between a supernova blast wave and inhomogeneous 

26 interstellar clouds formed by thermal instability generates multiple reflected shocks, which 

27 can further energize cosmic-ray particles originally accelerated at the blast-wave shock 
2B and produce the spectral break. Since the gamma-ray bright regions are expected to be 

29 different in the aforementioned models, studying SNRs with large apparent sizes can help 

30 to disentangle the origin of the spectral features. 

31 The Cygnus Loop (G74. 0—8.5) is one of the most famous and well-studied middle-aged 

32 SNRs. The size (~ 3°) makes it an ideal candidate for detailed morphological studies in 

33 high-energy gamma rays since it is larger than the LAT angular resolution above a few 

34 hundred MeV. The large angular offset from the Galactic plane {b ~ —8°. 5) reduces the 
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problems of background contamination and permits a detailed study of the environment 
around the shock region by means of infrared/optical/UV observations. The shell-like X-ray 
emission from thermal plasm a is prorninent in the northern region of the remnant, with a 



3B blowout in the southern rim ( 



shells is non-thermal ( iKeen et al 



Kuetal 



19841 ) ■ The radio spectrum from the limb-brightened 



19731 ) ■ No correlation with dense molecular clouds 



40 has been reported, although blast waves on the western limb might encounter molecular 



material (jScoville et al 



19771 ). The distance from the Earth is estimated to be 540^80° pc 



based o n the proper mot ion of optical filaments in conjunction with models of non-radiative 



shocks (IBlair et al 



20051 ) ■ The age h as been estimated t o be ~ 2 x 10'^ yr based on plasma 



parameters derived from X-ray data (IMivata et al. 



shock model and X-ray measurements (ILevenson et al. 



19941) an d ~ 1.4 x 10^ yr based on the 



19981). 



Four LAT sources positionally associated with the Cygnus Loop SNR are listed in the 



IFGL catalog (lAbdo et al. 



3) 



2010dl ). In this paper, we report a detailed analysis of Fermi 



48 LAT observations in the Cygnus Loop region. First, we briefly describe the observations 

49 and data selection in Section O The analysis procedure and the results are presented in 

50 Section [3l with the study of the morphology and spectrum of emission associated with the 

51 Cygnus Loop. Results are then discussed in Section H] and our conclusions are presented in 

52 Section [5l 



2. OBSERVATIONS AND DATA SELECTION 



54 The LAT is the main instrument on Fermi, detecting gamma rays from ~ 20 MeV 

55 to > 300 GeVlj. Details abou t the LAT inst r umen t and pre-launch expectations for the 



56 performance can be found in 



Atwood et al. 



(120091 ). Compared to earlier high-energy 



^ As noted below in the present analysis we use only events with energies > 200 MeV. 
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57 gamma-ray telescopes, the LAT has a larger field of view (~ 2.4 sr), a larger effective 

5B area (~ 8000 cm^ for >1 GeV on-axis) and improved point-spread function (PSF; the 68% 

59 containment angle > 1 GeV is smaller than 1°). 

60 Routine science operations with the LAT began on August 4, 2008. We have analyzed 

61 events in the region of the Cygnus Loop collected from August 4, 2008, to August 1, 2010, 

62 with a total exposure of ~ 6 x 10^° cm^ s (at 1 GeV). The LAT was operated in sky-survey 

63 mode for almost the entire period. In this observing mode the LAT scans the whole 

64 sky, obtaining complete sky coverage every 2 orbits (~ 3 hr) and approximately uniform 

65 exposure. 

66 We used the standard LAT analysis software, the ScienceTools version v9rl6, publicly 

67 available from the Fermi Science Support Center (FSScfl and applied the following event 

68 selection criteria: 1) events have the highest prob ability of being garn ma rays, i.e., they 



20091 ). 2) events have a 



are classified in the so-called Pass 6 Diffuse class (lAtwood et al. 
reconstructed zenith angle less than 105°, to minimize the contamination from Earth-limb 
gamma-ray emission, and 3) only time intervals when the center of the LAT field of 
view is within 52° of the local zenith are accepted to further reduce the contamination 
by Earth's atmospheric emission. We also eliminated two short periods of time dur ing 



74 which the LAT detected the bright G eV-emitting GRB 081024B flAbdo et al 



75 GRB 1001 16A f McEnerv et al 



2010d) and 



2010l ) within 15° of the Cygnus Loop. We restricted the 



analysis to the energy range > 200 MeV to avoid possible large systematics due to the 
rapidly varying effective area and much broader PSF at lower energies. 



^Software and documentation of the Fermi ScienceTools are distributed by Fermi Science 
Support Center at http://fermi.gsfc.nasa.gov/ssc 



- 7- 



, ANALYSIS AND RESULTS 
3.1. Morphological analysis 
3.1.1. Method 



In order to study the morphology of gamma-ray emission associated with the Cygnus 

a 



82 Loop we per: 



Mattox et al. 



ormed a binned likelihood analysis based on Poisson statisticc|(see e.g. 



19961 ). We used only events above 0.5 GeV (compared to the 0.2 GeV used in 

84 the spectral analysis) for the morphological study to take advantage of the narrower PSF 

85 at higher energies. For this work we used the instrument response functions (IRFs) P6_V3, 

86 which were developed following the launc h to address gamma -ray detection inefficiencies 



that are correlated with background rates (IRando et al 



20091 ) . The analysis was performed 



88 over a square region of 12°xl2° width with a pixel size of 0.°1. We set the centroid of 

89 the region to (R.A., Dec.) = (2l'^03'"03^ 33°42'56"), 2° shifted from that of the Cygnus 

90 Loop toward negative Galactic latitudes to avoid the background given by Galactic diffuse 

91 emission and Galactic sources. Figured] (a) shows a count map in the 0.5-10 GeV energy 

92 band in the region used for the analysis, as well as the position of the Cygnus Loop 

93 from radio measurements and point sources in the IFGL catalog. The four LAT sources, 

94 IFGL J2046.4+3041, IFGL J2049. 1+3142, IFGL J2055. 2+3144 and IFGL J2057.4+3057, 

95 are associated with the Cygnus Loop. Note that no gamma-ray pulsation was found for any 

96 of these LAT sources. 



■^As implemented in the publicly available Fermi Science Tools. The documenta- 
tion concerning the analysis tools and the likelihood fitting procedure is available from 



http: / /fermi. gsfc.nasa.gov/ssc/data/analysis/documentation/Cicerone/ 
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3.1.2. Background model 



9B Although the Cygnus Loop is at intermediate Galactic latitude, the contribution of 

99 the Galactic interstellar emission in the gamma-ray band is still important; it must be 

100 carefully modeled to perform morphological studies. Some of the interstellar gas tracers in 

101 the standard diffuse model provided by th e LAT collaboration are not fully adequate for the 



Cygnus region, notably the E(B-V) map (ISchlegel et al.l 119981 ) because of infrared source 



103 contamination and temperature correction problems in such a massive-star forming region. 

104 We therefore constructed a dedicated diffuse emission model. The model is analogous 

105 to the standard LAT diffuse model and it includes: a) an isotropic background, taking 

106 into account the isotropic diffuse gamma-ray emission as well as residual misclassified CR 

107 interactions in the LAT; b) large-scale Galactic inverse Compton emission produced by CR 
lOB electrons and positrons upsc attering low-energy photons, modeled using the GALPROP 



code (e.g. 



Porter et al. 



20081 ): c) emission from interstellar gas arising from nucleon-nucleon 
no interactions and electron Bremsstrahlung, which is modeled through spatial templates 

111 accounting for atomic gas and CO-bright molecular gas, partitioned along the line of sight 

112 to separate the Cygnus complex from the segments of the spiral arms in the outer Galaxy 

113 seen in this direction, as well as dark gas traced by visual extinction. With respect to th e 



114 standard diffuse model, this one includ es higher-resolution 



115 visual extinction as a dark-gas tracer ( iRowles fc Froebrich 



2009 



I data (iTavlor et al. 



20031) 



Froebrich &: Rowles 



20101) and it is specifically tuned to reproduce LAT data in the Cygnus region, including 



117 the region of the Cygnus Loop. All these components, along with individual sources, were 

118 jointly fitted to the LAT data in 10 energy bands over the range 0.1-100 GeV with a free 

119 normalization in each energy bin (except for the inverse Compton mo del that was kept 



120 fixed). For further details we refer the reader to the dedicated paper (lAckermann et al. 



201ll ). where the model is also discussed in detail in terms of CR and interstellar medium 
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122 properties. We note that the presence of the Cygnus Loop was taken into account in this 

123 study. Several models were considered for the Loop, a combination of point sources and 

124 geometric templates such as a disk and a ring, as in the analysis performed in this paper. 

125 In this way we verified that the impact of the emissi on from the C ygnus Loop on the 



126 parameters of the global model of the region is small (ITibaldo 



201 ih . 



127 The results of this analysis were used to construct two model cubes, as a function 

128 of direction and energy, separately accounting for the isotropic and smooth large-scale 

129 Galactic inverse-Compton emission (a and b) and the structured emission from the gas (c). 

130 Such model cubes are part of the background model used to study the Cygnus Loop in 

131 this paper. For each of them we included a free normalization in order to further allow the 

132 model to adapt in the different cases we investigated along the paper. 

133 In addition to interstellar emission, the background model to study the Cygnus Loop 

134 includes individual point-like sources in the IFGL catalog within 15° of the Cygnus Loop 

135 except for the sources associated with the Cygnus Loop itself in the catalog; their positions 

136 were kept fixed at those given in the catalog and the spectra of the two gamma-ray pulsars 

137 in the region used for the analysis were modeled as power laws with exponential cutoffs 
13B leaving all spectral parameters free, while the spectra of the other sources were modeled as 

139 power laws leaving the integral fluxes as free parameters and assuming the spectral indices 

140 reported in the catalog. Note that, due to the PSF, which is poor compared to other 

141 wavelengths and strongly energy dependent, and the presence of a bright and structured 

142 background given by the interstellar emission, it is difficult to mask the background 

143 sources, and they are instead modeled along with the Cygnus Loop. The resulting model 

144 of background emission (i.e. not including emission associated with the Cygnus Loop) is 

145 shown in Figure [1] (b). The pulsars J2043+2740 and J2055+25 are the most important 

146 point sources in the vicinity, but the amount of events from those sources that fall within 
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147 the Cygnus Loop (due to the broad low-energy PSF) is only 0.4 % and 0.2 % of the 

148 estimated emission from the Loop, respectively. 



149 3.1.3. Comparison with observations at other wavelengths 

150 Figure [2] shows the count map after subtracting the background emission in a 6° x 6° 

151 region centered on the Cygnus Loop, (R.A., Dec.) = (20'*51"^06^, 30°41'00"), with overlays 

152 of images at different wavelengths: X-rays, Ha line, radio continuum at 1420 MHz, infrared 

153 radiation at 100 fim and CO 2.6 mm line. The correlation between gamma rays. X-rays 

154 and Ha emission is evident. There is correlation among gamma rays and radio continuum 

155 emission in the northern part of the Cygnus Loop. On the other hand, the southern rim 

156 is brighter in radio continuum emission than in gamma rays, a phenomenon that perhaps 

157 might be explained b y the existence of an other SNR overlapping with the southern part 



158 of the Cygnus Loop (lUyaniker et al. 



2OO2I ). The CO line intensities were integrated for 



No 



159 velocities with respect to the local standard of rest —25 km s~^ < < 30 km s' 

160 obvious association with CO emission is found; some molecular material is apparently 

161 located on the Western side of the Cygnus Loop, but the relationship is not clear. On the 

162 other hand, some correlation with thermal emission from dust at 100 fim, which can be 

163 considered as a proxy of total interstellar matter densities, is possible. 

164 To quantitatively evaluate the correlation with emission at other wavebands, we fitted 

165 the LAT counts with the different models for the Cygnus Loop on top of the background 

166 model described above. First the Cygnus Loop was modeled with the four IFGL sources, 

167 and then using the images at other wavelengths as spatial templates assuming a simple 

168 power-law spectrum. Note that we did not use the CO and infrared images as spatial 



^The velocity corresponding to the distance from the Earth (540 pc) is ~ km s ^. 
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169 templates due to clear differences between them and the gamma-ray image as shown in 

170 Figure [H The resulting maximum likelihood values with respect to the maximum likelihood 

171 for the null hypothesis (no emission associated with the Cygnus Loop) are summarized in 



172 Table [H The test statistic (TS) values, i.e. — 2 ln(likelihood ratio) (e.g. iMattox et al.lll996l ). 

173 for the X-ray and Ha images are significantly larger than for the four IFGL sources. On 

174 the other hand, the TS for the radio image increases moderately in spite of the association 

175 in the northern rim, confirming that radio continuum structures in the southern rim do not 

176 well correlate with gamma-ray emission. 



3.1.4- Geometrical Models 



178 We further characterized the morphology of gamma-ray emission associated with 

179 the Loop by using simple parametrized geometrical models. We started with a uniform 

180 disk/ring assuming a simple power-law spectrum. We varied the radius and location of the 

181 disk and evaluated the maximum likelihood values. In the case of the ring, we varied inner 

182 and outer radii as well. The resulting TS values are reported in Table [H The TS value for 

183 the ring with respect to the disk shape is ^ 12. Assuming that, in the null hypothesis, 

184 the TS value is distributed as a with n degrees of freedom, where n is the difference in 

185 degrees of freedom between the two nested models compared □ (n = 1 in the present case), 

186 it would be equivalent to an improvement at ~ 3.5 a confidence level. Let us note, however, 
that the conversion of TS values into co nfidence level (or , equiy alently, false positive rate) 



is subject to numerous caveats, see e.g. 



Protassov et al. 



(I2OO2I). We will thus take into 



189 account the source morphology uncertainties in the spectral analysis, below. In order to 

190 further illustrate the morphology of the gamma-ray emission, in Figure [3] we show its radial 



^see link to Fermi Science Tools Cicerone 
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191 profile compared with the best-fit disk/ring models. 

192 Finally, we want to verify if there are any spectral variations in the gamma-ray emission 

193 associated with the Cygnus Loop we are modeling as a whole. We thus divided the best-fit 

194 ring into four regions as shown in Figure H] and allowed an independent normalization and 

195 spectral index for the four portions of the ring. There was no significant improvement of 

196 the likelihood for such a non-uniform ring. The TS value and power-law spectral index for 

197 each of the four regions of the remnant are reported for reference in Table [2l No significant 

198 differences are found between the four spectral indexes. Therefore we adopted the uniform 

199 ring template with maximum likelihood parameters for the whole SNR in the following 

200 spectral analysis. 



201 3.2. Spectral analysis 

202 To measure the spectrum we made maximum likelihood fits in 8 logarithmically-spaced 

203 energy bands from 0.2 GeV to 100 GeV, using the ring template as the model for the spatial 

204 distribution of the Cygnus Loop. Figure [5] shows the resulting spectral energy distribution 

205 (SED). Upper limits at the 90 % confidence level are calculated assuming a photon index of 

206 2 if the detection is not significant in an energy bin, i.e., the TS value with respect to the 

207 null hypothesis is less than 10 (corresponding to 3.2 cr for one additional degree of freedom). 

208 Note that the value of the spectral index has a negligible effect on the upper limits. 

209 We identify at least three different sources of systematic uncertainties affecting the 

210 estimate of the fluxes: uncertainties in the LAT event selection efficiency, the morphological 

211 template and the diffuse model adopted for analysis. Uncertainties in the LAT effective 

212 area are estimated to be 10 % at 100 MeV, decrea sing to 5 % at 500 MeV, and increasing 



213 to 20 % at 10 GeV and above (IRando et al. 



20091 ). Evaluating the systematic uncertainties 
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214 due to the modeling of interstellar emission is a challenging task, because interstellar 

215 emission is highly structured and methods used at other wavelengths, like comparisons with 

216 neighboring regions, are not fully adequate in the GeV band. We therefore roughly gauged 

217 the related uncertainties by comparing the results with those obtained by adopting instead 

218 the standard LAT diffuse background models 0. We similarly gauged the uncertainties due 

219 to the morphological template by comparing the results with those obtained by using the 

220 best-fit disk template instead of the ring. The total systematic errors are set by adding the 

221 above uncertainties in quadrature. Systematic uncertainties are driven by the imperfect 

222 knowledge of the background emission and, especially below a few hundred MeV, of the 

223 LAT response. In Figure [5] we show the uncertainties obtained following these prescriptions. 

224 We probed for a spectral break in the LAT energy band by comparing the likelihood 

225 values of a spectral fit over the whole energy range considered based on a simple power 

226 law and other spectral functions. Note that no systematic uncertainties are accounted for 

227 in the likelihood fitting process. The TS values and best-fit parameters are summarized 

228 in Table [3l The fit with a log-parabola function yields a TS value of ~ 50 compared to a 

229 simple power-law model, which corresponds to an improvement at the ~ 7 cr confidence 

230 level. In spite of the uncertainties discussed above in the estimate of the confidence level, 

231 the large TS value is indicative of a significant improvement in the fit. A smoothly broken 

232 power law provides a very slight increase in the likelihood with respect to the log-parabola 

233 function, while a power law with exponential cutoff gives a worse fit. In conclusion, a 

234 simple power law as spectral model can be significantly rejected and from all the different 

235 models with cutoffs we get evidence for a steepening of the spectrum above 2-3 GeV. 

236 We detect gamma-ray emission with a formal significance of 23 a for the above curved 

^ gll_iem_v02 and isotropic-iem-v02 available from the FSSC 
http : / / f ermi . gsf c . nasa . gov /ssc/ dat a/access/lat /B ackgroundMo dels . html 
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237 spectral shapes. The observed photon flux and energy flux in the 0.2-100 GeV range are 

238 5.0lo;g X 10"^ cm"^ s~^ and x 10"^^ erg cm~^ s~^, respectively. 



DISCUSSION 



240 The gamma-ray luminosity inferred from our analysis is ~ 1 x lO^^erg s ^ between 

241 1-100 GeV, lower by one order of magnitude than obse rved for other GeV-emitting 



242 SNRs (typically > 10^^ erg s'^, 



Abdo et al 



2009 



2010a 



J). The spatial distribution is 



243 best represented by a ring with inner/outer radii 0°.7 ± 0°.l and 1°.6 ± 0°.l, respectively. 

244 This makes the Cygnus Loop the largest gamma-ray emitting SNR observed so far, allowing 

245 us to perform a detailed morphological comparison with emission at other wavelengths. 

246 There is a correspondence among gamma-ray emission. X-ray rims and Ha fllaments, 

247 indicating that the high-energy particles responsible for gamma-ray emission are in the 
24B vicinity of the shock regions. The Balmer-dominated filaments define the current location 

249 of the blast wave and mark the presence of neutral material. Detailed studies of the 

250 particul ar locations at the northeast haye used these nonradiatiy e shocks as density 



251 probes ( iRaymond et al. 



1983 



Long et al. 



1992 



Hester et al 



1994} and derived post-shock 



252 densities of ~ 5 cm ^ where gamma-ray emission is expected to be bright due to the 

253 compressed material and high density of accelerated particles. 

254 The radio continuum emission, originated by high-energy electrons via synchrotron 

255 radiation, is well correlated with gamma-ray emission in the northern region of the remnant 



256 but no t in the southern one The presence of a second SNR was suggested by 



Uyaniker et al. 



257 ( 120021 ). The two SNRs wo uld be at about t he same distance based on the rotation measure 



25B analysis of the radio data (jSun et al. 



20061 ). The lack of correlation between gamma rays 



259 and radio continuum emission in the southern region plausibly implies that the second SNR 
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260 is not producing significant gamma-ray emission at our current sensitivity. There might be 

261 some correlation between total matter densities as traced by infrared thermal emission from 

262 dust and gamma-ray emission, whereas CO emission does not obviously overlap with the 

263 Cygnus Loop. 

264 From these considerations, we argue that the bulk of gamma-ray emission comes from 

265 interactions of high-energy particles accelerated at the shocks of the Cygnus Loop with 

266 interstellar matter or fields in the regions just adjacent to the shocks with a gas density of 

267 ~ 5 cm 

268 To model the broadband emission from the entire SNR we adopt the simplest possible 

269 assumption that gamma rays are emitted by a population of accelerated protons and 

270 electrons distributed in the same region and characterized by constant matter density and 

271 magnetic field strength. We assume the injected electrons to have the same momentum 

272 distribution as protons. This assumption requires a break in the momentum spectrum 

273 because the spectral index in the radio domain, corresponding to lower particle momenta, 

274 is much harder than for gamma rays, which correspond to higher particle momenta. 

275 Therefore, we use the following functional form to model the momentum distribution of 

276 injected particles: 

^ 9"! -{sh-sl)/2 

277 where pbr is the break momentum, sl is the spectral index below the break and sh above 

278 the break. Note that here we consider minimum momenta of 100 MeV since the details 

279 of the proton/electron injection process are poorly known. 

280 Electrons suffer energy losses due to ionization. Coulomb scattering, Bremsstrahlung, 

281 synchrotron emission and inverse Compton (IC) scattering. We calculated the evolution of 
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282 the electron momenta spectrum by the following equation: 



e,p 



dt 



d_ 

dp 



(2) 



where h^^p = —dp/dt is the momentum loss rate, and is the particle injection rate. 
We assume Q^p to be constant, i.e., that the shock produces a constant number of 
particles until the SNR enters the radiative phase, at which time the source turns off. This 
prescription approximates the weakening of the shock and the reduction in the particle 



287 acceleration efficier icy, which would 



compression ratio (IMoraal fc Axford 



3e pro perly treated by using a time-dependent shock 



19831 ). To derive the remnant emission spectrum we 
calculated Ne^p{p,To) numerically, where Tq is the SNR age of 2 x 10'^ yr. Note that we 
neglected the momentum losses for protons since the timescale of neutral pion production 
is ~ 10''/nH yr where nn is the gas density averaged over the entire SNR shell and is 
much longer than the SNR age. Also we do not consider the gamma-ray emission by 
secondary positrons and electrons from charged pion decay, because the emission from 
secondaries is generally unimportant relative to that from primary electrons unless the gas 
density is as high as that in dense molecular clouds and the SNR evolution reaches the 
later stages, or the injected electron-to-proton ratio is much lower than locally observed. 
The gamma-ray spectrum from 7r° decay produced by the interactions of protons with 



298 ambient hydrogen is calculated based on 



Dermerl (119861 ) using a scaling fa ctor o. 



299 account for helium and heavier nuclei in target material and cosmic rays (IMori 



1.84 to 



20091). 



300 Contributions from bre msstrahlung and inver s e Com pton scattering by accelerated electrons 



ar e computed based on 



302 on 



B 



Crusius fc Schlickeisei 



umenthal fc Gould! (Il970l ). whereas synchrotron radiation is based 

(igsj 



303 First, we consider a vr^-decay dominated model. The number index of protons in the 

304 high-energy regime is constrained to be sh ~ 2.6 from the gamma-ray spectral slope. The 

305 spectral index of the proton momentum below the break is determined to be sl ~ 1.8 by 
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306 modeling the radio spectrum as synchrotron radiation by relativistic electrons (under the 

307 assumption that protons and electrons have i dentical injection spe ctra). The spectral index 



20041 ). where a is defined as 



308 a of the radio continuum emission is ~ 0.4 (lUyaniker et al. 

309 Sy oc with Sy and v the fiux density and the frequency, respectively. It is difficult to 

310 derive the break momentum of the proton spectrum from the gamma-ray spectrum, since in 

311 the GeV energy band we expect a curvature due to kinematics of 7r° production and decay. 

312 The gamma-ray spectrum provides thus only an upper bound for the momentum break at 

313 ~ 10 GeV c~^. On the other hand, the momentum break cannot be lower than ~ 1 GeV 

314 to avoid confiicts with radio data. We adopt a break at the best-fit value, 2 GeV c~^. The 

315 resulting total proton energy, ~ 2.6 x lO''^ • (5 cm^^/nn) • (c?/540 pc)^ erg, is less than 

316 1 % of the typical kinetic energy of a supernova explosion. For an electron-to-proton ratio 

317 i^ep = 0.01 at 1 GeV c~^, which is the ratio measured at the Earth, the magnetic field 
31B strength is constrained to be -B ~ 60 |J.G by radio data. The magnetic field st rength of the 



319 undisturbed medium in the northeastern rim was estimated to be ~ 20 \xG by 



van der Laan 



320 (119621 ) based on the measurements of shell thickness and expansion velocities together 



321 with the theory of hydro magnetic shock pro pagation given the density of the undisturbed 



322 medium ~ 1 cm 



le-g. 



Hester et al. 



1994j ). The compression behind the shock front 



323 indicates a magnetic field strength similar to the value used above in the modeling. Using 

324 the parameters summarized in Table HI we obtained the SEDs shown in Figure [6] (a). 

325 It is difficult to model the gamma-ray spectrum with a model dominated by electron 

326 bremsstrahlung because the break in the electron spectrum required to reproduce the 

327 gamma-ray spectrum would appear in the radio domain as shown in Figure [6] (b). 

328 The gamma-ray spectrum can be reproduced by an inverse Compton dominated 

329 model shown in Figur^E] (c). Gamma-ray emission of IC origin is due to interactions of 

330 high-energy electrons with optical and infrared radiation fields and the cosmic microwave 
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background (CMB). We used in our calculations the first two components as they are 



modeled in 



Porter et al 



( 120081 ) at the location of the Cygnus Loop. Since their spectra 
are very complex, they are approximated by two infrared and two optical blackbody 
components. The flux ratio between the IC and the synchrotron components constrains the 
magnetic field to be less than 2 |j.G and requires a low gas density of tt-h ~ 2 x 10~^ cm~^ 
to suppress the electron bremsstrahlung. Alth ough such a low density may exist inside 



Ku et al. 



1984 ). gamma-ray emission peaks 



337 the remnant based on X-ray observations (e.g. 

33B at the shock regions where the gas density is ~ 1 — 5 cm ~^ (see above). Increasing the 

339 intensity of the interstellar radiation field would loosen the constraint on the gas density. 

340 However, a radiation field about 50 times more intense is required to satisfy the above 

341 assumption on the gas density. 

342 To summarize, it is most natural to assume that gamma-ray emission from the Cygnus 

343 Loop is dominated by decay of 7r° produced in nucleon-nucleon interactions of hadronic 

344 cosmic rays with interstellar matter. It should be emphasized that our observations of 

345 the Cygnus Loop combined with the radio data constrain the proton momentum break 

346 to be in the range, 1-10 GeV c~^, despite the lack of association with dense molecular 

347 clouds unlike the other middle-aged SNRs detected with the LAT. Thus in this case 
34B cosmic rays responsible for gamma-ray emission are localized near their acceleration 

349 sites without significant diffusion taking place. The correspondence observed between 

350 gamma rays and Ha e mission may be accounted for in the "crushed cloud" scenario by 



Uchiyama et al. 



( l2010l ). although the expe cted filamen t s can not be resolved by current 



352 gamma-ray telescopes. The predictions by 



Inoue et al 



(120 lOf ) cannot be directly compared 



353 to the Cygnus Loop since their simulations were performed for environments characterized 

354 by dense clouds. However, the scenario of acceleration by reflect ed shocks might be 



355 operative, on consideration of X-ray and optical observations (e.g. 



Graham et al. 



19951). 



- 19 - 



5. CONCLUSIONS 

357 We analyzed gamma-ray measurements by the LAT in the region of the Cygnus Loop, 

35B detecting significant gamma-ray emission associated with the remnant. The gamma-ray 

359 luminosity is ~ 1 x 10^'^erg s~^ between 1-100 GeV, lower than for other GeV-emitting 

360 SNRs studied with LAT data. The morphology of gamma-ray emission is best represented 

361 by a ring with inner/outer radii 0.°7 ± 0°.l and l.°6 ± 0°.l. The Cygnus Loop is 

362 thus the most extended gamma-ray emitting SNR detected in the GeV band so far and 

363 the morphology of gamma-ray emission can be compared in detail with observations at 

364 other wavelengths. There is correspondence among gamma rays, the X-ray rims and the 

365 Ha filaments, indicating that the high-energy particles responsible for the gamma-ray 

366 emission are in the vicinity of the shock regions. 

367 The gamma-ray spectrum has a break in the 2-3 GeV energy range. The decay 
36B of 7r° produced by interactions of hadrons accelerated by the remnant with interstellar 

369 gas naturally explains the gamma-ray spectrum. In this scenario our observations of the 

370 Cygnus Loop indicate that the proton momentum spectrum is steep in the high-energy 

371 regime, with a spectral break which is constrained together with radio continuum emission 

372 in the range 1-10 GeV c~^. The absence of molecular clouds in the areas of gamma-ray 

373 emission (contrary to other middle-aged Fermi SNRs) constrains some of the scenarios 

374 invoked to explain the observed spectral properties of GeV emitting SNRs. 
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This manuscript was prepared with the A AS lAT^X macros v5.2. 



Galactic longitude 



6 8 

[counts/pixel] 



Fig. 1. — (a) Fermi LAT count map in the Cygnus Loop region for photon energies 0.5- 
10 GeV. The count map has a pixel size of 0.°1 and is smoothed for display with a Gaussian 
kernel of a = 0.°5. Note that all along the paper the analysis is conducted on unsmoothed 
data taking into account the instrument PSF in the likelihood analysis. The white circle is 
the location of the Cygnus Loop, defined by its radio emissi on. Green crosses indicate the 



2010dh . (b) count 



positions of gamma-ray sources listed in the IFGL catalog ( lAbdo et al. 
map expected from the background model (taking into account the LAT PSF). The four 
LAT point sources associated with the Cygnus Loop are not included in the model. The 
image is binned and smoothed in the same manner as the real data. 
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Fig. 2. — Background-subtracted LAT count map in the 0.5-10 GeV energy range. The count map is 
binned using a grid of 0.°05 and smoothed with a Gaussian kernel of cr = 0.°5. Negative residuals are shown 
to gauge the quality of the subtraction of the background emission. Green contours correspond to images at 
different wavelengths, (a) X-ray count map (0.1-2 keV) by ROSAT. Contours are at 20, 40, 60, 80 % levels; 
the image was first cleaned from background emission, estimated by fitting data surrounding the Cygnus 
Loop with a bilinear function, and smoothed with a Gaussian kernel of cr = 0.°2; (b) Ha image obtained from 
the publicly available Digital Sky Survey obtained with the same procedure explained for X-ray data. We 
selected the POSS-II F (r ed) filtered s urvey whose transmission coefficient peaked near Ha . (c) 1420 MHz 



radio continuum emission ( Reich 



Il982[ ): extraction of the contours as for the previous images, (d) ^-^CO (J 
1 — > 0) line inte nsities integrated f or velocities from —25 km s^^ to 30 km s The data are take n from 



the C fA survey (jPame et al 



20011) cleaned from background using the moment-masking technique (jPame 



20111 ): the image was smoothed using a Gaussian kernel with a = 0.°25; contours are at 1,4, 7, 10 K km s ^ 



(e) The infrared intensity map at 100 fim by InfraRed Astronomical Satellite (IRAS) (jBeichman et al.ll 19881 ): 



the image was smoothed using a Gaussian kernel of <t = 0.°2; contours are at 15, 25, 35, 45 MJy sr^^. The 
contour at the top-right corner is the highest one. (f) the effective LAT PSF in the energy band of the LAT 
count map for a photon spectral index of 2.5. The PSF map is binned and smoothed in the same manner as 
the real data. 



0.5 1 1.5 2 2.5 

eO 

Radial profile of the Cygnus Loop in gamma rays in the 0.5-10 GeV energy 
range (crosses). The origin is the center of the best-fit ring model. Gamma-ray data have the 
background emission subtracted. Note that the data are not smoothed. Overlaid are the dis- 
tributions expected for the best-fit ring shape (solid line) and the best-fit disk shape (dotted 
line) as emission surfaces, with parameters fit to gamma-ray data taking the LAT instrument 
response into account. Details of the fits are described in the text. 



3 



Fig. 3. 
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Fig. 4. — Definitions of tlie regions of the Cygnus Loop used for the morphology analysis 
(§ I3.ip overlaid on the LAT count map as shown in Figure [2J 
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Fig. 5. — Spectral energy distribution of the gamma-ray emission measured by the LAT for 
the Cygnus Loop. Red squares are LAT flux points. Horizontal bars indicate the energy range 
the flux refers to. Vertical bars show statistical errors in red and systematic errors (added 
in quadrature for illustration purposes) in black. In energy bins where the detection is not 
significant (test statistic < 10) we show upper limits at the 90 % confidence level. The blue 
region is the 68 % confidence range (no systematic error) of the LAT spectrum assuming 
that the spectral shape is a log parabola. 
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Fig. 6. — Multi-band spectrum of the Cygnus Loop. In t he GeV band LAT ni easurements 



20041 ) is shown by 



are reported as in Figure [51 The radio continuum emission (jUyaniker et al.l 
black dots. Radio emission is modeled as synchrotron radiation, while gamma-ray emission is 
modeled by different combinations of vr'^-decay (long-dashed curve), bremsstrahlung (dashed 
curve), and inverse Compton (IC) scattering (dotted curve). Details of the models are 
described in the text: a) Tr'^-decay dominated model, b) bremsstrahlung dominated model, 
c) IC dominated model. 
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Table 1: Test Statistics for Different Spatial Models Compared with the Null Hypothesis of 
No Gamma-ray Emission Associated with the Cygnus Loop (0.5-100 GeV) 



Model 


Test Statistic'^ 


Additional Degrees of Freedom 


Null hypothesis'^ 








4 point sources'^ 


318 


8 


/ZO^^T X-rays (0.1-2keV)'i 


406 


2 


Ha ^ 


434 


2 


1420MHz radio continuum*^ 


343 


2 


Uniform disk^ 


441 


5 


Uniform ring^ 


453 


6 


Non-uniform ring^ 


464 


12 



^—2 ln{Lo/L), where L and Lq are the maximum likehhoods for the model with/without the somce component, 
respectively. 

''Background only (no model for the Cygnus Loop). 



2010dD 



■^The four sources listed in the IFGL source list associated with the Cygnus Loop (|Abdo et al. 
'^Background-subtracted as described in Figure [21 

''The best fit parameters are: radius 1°.7 ± 0°.l and centroid (R.A., Dec.) = (20''52™, 30''50'). The error of 
the centroid is 0°.04 at 68 % confidence level. 

^The best fit parameters are: inner/outer radii O'^.T ± O'^.l, and 1°.6 ± 0°.l, centroid (R.A., Dec.) = (20''51™, 

30''50')- The error of the centroid is 0''.04 at 68 % confidence level. 
^The best-fit ring were divided into four regions as shown in Figure|3]and allowed an independent normalization 

and spectral index for the four portions of the ring. 
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Table 2: Test Statistics and Power-law Spectral Indexes for the Four Regions of the Remnant 

as Defined in Figure 1 ( 0.5-100 GeV) 

Region Test Statistic^ Spectral Index 

I 143 2.49 ± 0.10 

II 73 2.32 ± 0.12 

III 64 2.25 ± 0.15 

IV 41 2.37 ± 0.14 



2 ln{Lo/L), where L and Lq arc the maximum likchhoods for the model with/without the source component, 
respectively. 
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Table 3: Test Statistics and Parameters for Various Spectral Models (0.2-100 GeV) 



Spectral Model 


Test Statistic'^ 


Degrees 
of Freedom 


Spectral Parameters 


Power law 





2 


E-P; p= 2.23 ±0.02 


Power law with 


42 


3 


E-^exp (-f); 


exponential cutoff 






p= 1.57±0.12 
= 3.02 ± 0.65 GeV 


Log Parabola 


50 


3 


/ E \ PI P2log(, g^v) 
\1 GeV/ 

pi = 2.02 ± 0.03 
P2 = 0.27 ±0.02 


Smoothly broken power law 


51 


4 





Pi = 1.83 ±0.06 
P2 = 3.23 ±0.19 
= 2.39 ± 0.26 GeV 



^—2ln{Lo/L), where L and Lo arc the maximum likelihood values for the model under consideration and the 
power-law model, respectively. 

Note. — The test statistics for the best-fit uniform ring with exponential cutoff, log parabola, and smoothly 
broken power law with respect to the null hypothesis of no emission associated with the Cygnus Loop are 
572, 580, and 581 in the energy band 0.2-100 GeV. 
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Table 4: Model parameters for the Cygnus Loop. 



Model 
























(GeV c-i) 






(cm^'^) 


(10^8 erg) 


(1048 erg) 


(a) Pion 


0.01 


1.8 


2 


2.6 


60 


5 


2.6 


4.9 X 10~2 


(b) Bremsstrahlung 


1 


1.8 


2 


2.7 


12 


5 


0.21 


0.43 


(c) Inverse ComptonS 


1 


1.8 


25 


5.0 


1.8 


0.02 


5.9 


9.8 



'^Thc ratio electrons-to-protons at 1 GeV c ^ . 

'^The momentum distribution of particles is assumed to be a smoothly broken power-law, where the indices 
and the break momentum are identical for both accelerated protons and electrons, sl is the spectral index 
in momentum below the break. 

'^Phr is the break momentum. 

'^Spectral index in momentum above the break. 

"Average hydrogen number density of ambient medium. 



20051 ) ■ The total energy is calculated for 



^The distance from the Earth is assumed to be 540 pc (IBlair et al 
particles > 100 MeV cr^. 

^Seed photons for inverse Compton scattering of electrons include the CMB, two infrared (Tir ~ 34,4.7 x 
10^ K, Um = 0.34,6.3 x 10"^ eV cm~^, respectively), and two optical components (Topt ~ 3.6 x 10^,9.9 x 
10"^ K, Uopt = 0.45,0.16 eV cm~^, respectively) in the vicinity of the Cygnus Loop. 



